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Abstract 

We study the quantum dynamics of strongly interacting few-boson mixtures in one-dimensional traps. 
If one species is strongly localized compared to the other (e.g., much heavier), it can serve as an effec- 
tive potential barrier for that mobile component. Near the limit of infinite localization, we map this to a 
system of identical bosons in a double well. For realistic localization, the backaction of the light species 
on the "barrier" atoms is explained — to lowest order — in terms of an induced attraction between these. 
Even in equilibrium, this may outweigh the bare intra-species interaction, leading to unexpected correlated 
states. Remarkably, the backaction drastically affects the inter-species dynamics, such as the tunneling of 
an attractively bound pair of fermionized atoms. 
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Tunneling of particles through an energetically forbidden region is a hallmark quantum- 
mechanical effect, which illustrates the wave nature of matter. The experimental flexibility of 
ultracold atoms makes it possible to study this phenomenon in a clean and highly controlled en- 
vironment, where the role of the tunnel barrier is played, e.g., by light forces such as in optical 
lattices jhj] . This has fostered the direct observation of fundamental effects like second-order tun- 
neling |2], Josephson oscillations and nonlinear self-trapping yfl. Moreover, optical lattices have 
proven powerful quantum simulators, giving insight into the role of tunneling in, e.g., the quantum 
phase transition from superfluid to insulator [4], or in spin-exchange processes responsible for 
quantum magnetism [5I]. 

All of these cases reflect the paradigm of tunneling through a classical barrier described by 
some external potential. In this Letter, we investigate tunneling through a material or quantum 
barrier, in the sense that it interacts with the particles. This is realized via quasi-one-dimensional 
(ID) mixtures of two atomic species, one of which is squeezed in the trap center. When the lo- 
calization is very tight, we show that this can indeed be understood as an effective tunnel barrier 
for the other species. As that confinement is relaxed, the "barrier" atoms move due to the backac- 
tion of the other species. The dramatic effect of this correlation on the mobile species is studied 
both for the ground state and for the inter-species tunneling dynamics. To first explore its mi- 
croscopic mechanism, we consider a few-atom system system, which is studied in a numerically 
exact fashion utilizing the multi-configuration time-dependent Hartree method [6]. This wave- 
packet dynamics tool which has been applied successfully to systems of identical bosons as well 
as mixtures (see QV[] for details). We then go on to derive approximate models that capture the rele- 
vant physics for higher atom numbers. Still, such a small system may be achieved experimentally, 
e.g., by creating arrays of ID optical lattices (such as in Jst) or of ID tubes [sj], each contain- 
ing only a few atoms. Moreover, high-resolution imaging techniques such as scanning-electron 
microscopy allow for single-site addressability |^]. 

Model. — We consider a mixture of two ID bosonic species, labeled o = A, 5. These may 
correspond to different atomic species (or isotopes); however, in the case of equal masses they 
can also be thought of as different hyperfine components. The many-body Hamiltonian then 
reads H = £ CT H G + Hab , with the single-species Hamiltonian H = J\_f] Pa i + (xo,i) + 
Hi<j8aS(x ,i-x a j) and the inter-species coupling H AB = Lj^Ab5(x A) ( -x Bj ). In what fol- 
lows, we will focus on the case of harmonic trapping potentials U a (x) = ^m a co^x 2 and repulsive 
forces, gcr , £ab > 0. (Note that taking g a — > °°, the component o can be mapped to a fermionic one 
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llOll . referred to asfermionization. In this sense, considering bosonic atoms poses no restriction in 
ID.) Rescaling to harmonic-oscillator units, we can eliminate fi = niA = C0a = 1 fcn\. 

Material tunneling barrier. — Imagine the situation where one of the species (say B) is much 
more strongly localized in the trap center. This can be achieved by drastically reducing its (un- 
perturbed) length scale a# = 1 / \Jm B C0q <C 1, which amounts to having a near- zero mass ratio 
a = /ma/^b 1 (for different atom masses such as in Li/Cs) and/or strong confinement of B, 



CO a/ C0r <C 1 (e.g., realized via species-dependent optical lattices Ill2ll ). For simplicity, let us as- 
sume much heavier B atoms but equal frequencies. This is not crucial, though, and experimentally 
a large frequency cob may enhance the localization effect. 

Being tightly localized, the B atoms should feel no density variations of A and, to lowest order, 
the total density matrix « /3^ a) <g> p^ B) can be approximately factorized. (This even holds for 
strong repulsion gAB ^> 1 , since the strong localization of the B atoms at the trap center effectively 



reduces the interaction to a single-particle potential II 1411 .) Integrating out the heavy B atoms leads 
to an effective Hamiltonian for the light species [Q], 

=H K +tT B [H A Bp^ B) ] =H A + gAB^MxA,i). (1) 

i 

In this light, A ought to feel only an effective single-particle potential §ab^b {x), which in our case 
likens a sharp barrier at x = 0, with a width ~ a# given by the one-body density pB = n^/N^ (in 
terms of the number density n%) and its height proportional to the inter-species coupling gAB- In 
particular, for N#g ab 3> 1> this barrier practically splits the harmonic trap for the A atoms into a 
pronounced double well. (By extension, localizing the B atoms over several centers should make 
for an effective "lattice" for A.) This scenario naturally brings up the question whether one can 
create an effective localization of the light atoms (in equilibrium) or even see dynamic effects such 
as tunneling solely due to inter- species interactions. 

Figure \T\ shows that, near the limit of strong localization, the simple picture invoked above is 
indeed consistent with our exact results. In detail, we have evidenced this on the ground state 
of Na = 4 atoms in the presence of one localized "impurity" (a = 0.001). Absent any intra- 
species interactions (#a = 0), the A atoms are coherently spread over the left- and right-hand 



side of the trap in analogy to an external double- well potential II 1111 . even though the system as a 
whole is strongly interacting (gAB = 25 unless otherwise stated). This materializes in the two-body 
density Pa\ x \^ x i)^ measuring the joint probability density of finding two A atoms at x\ and X2: 
In Fig. [Da), both are equally likely to be on the same site (peaks on the diagonal, x\ « jc 2 ) as on 
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Figure 1: (color online) Strong-localization limit, a = m^/m^ = 0.001, for N A = 4 atoms in the presence 

(2) 

of a B atom localized at the center. Top: Two-body density p A {x\,X2) for gA = 0, 0.5, and 25 from left to 
right. Bottom: Density profile p A (x) for different gA- (All quantities in harmonic-oscillator units of A, see 
text.) 

opposite ones (x\ ~ — X2). As increases, the A atoms assume a "Mott-insulator"-type state with 
all A atoms more or less localized in either well [g A = 0.5; Fig.QJb)], and eventually fermionize 
[g A = 25; Fig. EDc)]. The checkerboard pattern emerging in the latter case can be understood 
simply as a gas of A-component bosons with hard-core repulsion (or, equivalently, noninteracting 
fermions) immersed in a double-well trap Jll, 13]: In each well, measuring one A-boson at x\ pins 
down the position of the remaining bosons to A^a — 1 discrete spots. This reflects in the density 
profile Pa(x), where A^a marked density maxima form as gA — ^ 00 (Fig-EQ bottom). 
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Figure 2: (color online) Beyond the static -barrier picture: Two-body density p A ~ {x\,X2) for {top) 
gA = 0, with increasing mass ratios a = 0.001,0.002,0.02 (from left to right); (bottom) gA = 0.5, 
a = 0.001, 0.12, 0.2. Same atom numbers as in Fig.Q] 

Beyond the static limit. — We have so far proceeded on the assumption that the B atoms 
were heavy enough to be frozen out completely, that is, treated as a classical potential. For less 
restrictive mass ratios, however, these are expected to move due to the backaction of the light 
atoms. The impact of this correlation becomes palpable in Fig. [2l where the ground-state evolution 
of '(xi,X2) is displayed as the mass ratio a is increased. For ,£4 = (Fig.O top), the previously 
uncorrected pattern develops into a strongly localized one (a = 0.02), where all A atoms are 
found exclusively on the same site; i.e., they cluster. By contrast, the anti-correlated "insulating" 
ground state observed at gA = 0.5 turns into a seemingly uncorrelated one as the B atom becomes 
less heavy (bottom). 

To understand this wealth of phenomena, it would be desirable to extend the effective Hamil- 
tonian (OQ) to higher orders in a > 0. This is nontrivial since H(a) in itself does not suggest 
a straightforward power series in a. Our basic procedure is sketched in what follows (details 
will be given elsewhere Iil4ll ). Given the strong localization of the B atom (A^ = 1 without 
loss of generality), we introduce center-of-mass and Jacobian coordinates relative to xb ~ 0. 
The resulting Hamiltonian affords an expansion in powers of a, which we can exploit so as 
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to derive an effective Hamiltonian for the A-atoms by tracing out out B, which leaves us with 
H A (a) = ftp + agAB^Lfc// [(xA,k-XA,i)n' B (x A> i)-n B (x Aj i)] + 0(a 2 ). The first term recovers 
the initial infinite-mass approximation; the first-order correction in turn may be understood as an 
additional external potential 8Ua (x) = —OCgAB^T^ [ m ' B (x) + n B (x)], plus an induced (nonlocal) 
interaction between two A atoms: 8Va [x\ , X2) = OCgAB [x\n' B (x2) + X2n' B (x\)] /2, which adds to the 
local one V A (x\,X2) =g A S(xi —X2). 

Given the effective multi-well (here: two sites) geometry experienced by the A atoms, it is 
tempting to explain some qualitative features from the perspective of simplified lattice models. 
Since, for a — > 0, the single-particle physics of species A is governed by h A = \p 2 + \x 2 + 
g A B n B(x), it is natural to use its Bloch-type eigenstates (which for a double well are delocalized 
anti-/symmetric functions that come in bands of doublets) as a basis for the many-body Hamil- 



tonian. In the spirit of the (bosonic) Hubbard model Ill5h . it is actually preferable to introduce 
Wannier-type functions localized on the left (right) site s in band /3 = 0, 1, . . . Expanding 
the many-body Hamiltonian H A {a) in terms of these Wannier functions and — in the case of suf- 
ficiently small interaction energies and deep wells — retaining only lowest-band and on-site terms 
yields an effective Bose-Hubbard model 

4 BH) = -J L *fa + - 1)- (2) 

Here —J = (w^\h A + 8U A (x)\wty) is the renormalized tunnel coupling of the lowest band, and 
u = (wf^® 2 \VA + 8Va\ws°^ 2 } = uq + 8u denotes the on-site interaction. While the correction to 
J tends to be small, the on-site renormalization due to 8Va can have a huge effect if it is on the 
same order as the a = term. A closer analysis Jl^j reveals that 8u < always, signifying an 
attractive induced interaction. There is an intuitive way of picturing this self-interaction: If an A 
atom sits in one well, then it will repel the barrier; thus the well becomes more spacious, making 
it energetically favorable to accommodate yet another A atom. 

Taking the induced attraction at face value, this casts a light on our results above. For gA = 
(Fig. |2l top), u = 8u < for a finite mass of B. This makes it plausible that the A atoms tend 
to cluster on same sites at a ~ 0.02, although being really noninteracting. Likewise, for small 
repulsion as in Fig. [2] (bottom; gA = 0.5), both terms may even cancel, uq + 8u ~ - this helps 
explain the seemingly uncorrelated pattern as a is increased. 

Although giving the right trend for weak interactions between the A atoms, the validity of © is 
actually much more limited than that of the effective Hamiltonian HpJ^a). In particular, the Bose- 
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Figure 3: (color online) Inter-species tunneling dynamics of Na = 2 bosons through a B atom: Relative A 
population of the right-hand side over time, ps.(t), for different mass ratios, a = 0.001 ( — ), a = 0.01 (- - 
-), and a = 0.12 (• • • ). (a) gA = 0, (b) gA = 0.5, (c) gA = 25; all gAB = 8- Insets: Probability piit) of finding 
two A atoms on the same side. 

Hubbard model breaks down for strong interactions comparable to the band gap. In that regime, 
rich multi-band effects can be found for stronger intra-species correlations, which we will discuss 
below in the context of the atoms' quantum dynamics. 

Inter-species tunnel dynamics. — We have so far investigated the equilibrium situation of 
species A in the presence of an effective "barrier" composed of a second, localized component. 
It would be thrilling to learn how this affects the quantum dynamics - for instance, can the light 
atoms tunnel through the heavy ones, and how does the barrier atoms' motion influence this inter- 
species tunnel dynamics? 
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To answer this question, let us investigate the time evolution of the A bosons loaded initially 
in, say, the left-hand side of the trap, with the barrier atoms (B) tightly centered. This can be done, 
e.g., by displacing the trap center of Ua(x) or by blocking the other half with a laser beam. Upon 
release, the A atoms may tunnel through species B, this way moving the barrier atoms, which in 
turn modifies the effective potential. To monitor the dynamics, we have recorded the percentage 
of A-atoms on the right, PR(t) = Jq pA(x;t)dx, in dependence of a. 

Figure H a) displays this population dynamics for Na = 2 noninteracting A bosons (gA = 0). 
For a static double well, we simply expect a Rabi-type oscillation of the population between left 
and right Jl6]. Indeed, for a = 0.001, this is what we find. Increasing the mass ratio, the clear sine 
mode gives way to a more complex, two-mode oscillation (a = 0.01), until for much larger values 
(a = 0.12) the tunneling slows down drastically, with only a tiny faster modulation on top. This is 



reminiscent of second-order tunneling well known from repulsive atom pairs |2,[l6l]. In fact, the 
effective Bose-Hubbard model © suggests that this corresponds to attractively bound pairs in a 
double well, where single-atom Rabi tunneling is highly suppressed ll 1711 . This line of reasoning is 



(2) 

supported by the pair (or same-site) probability pi{t) = Js XVX2 >q\ P a [x\,X2',t)dx\dx2, measuring 
how likely it is to find two A-atoms on the same site Mlal : This is plotted in the inset of Fig. Ha), 
indicating that the atom pair becomes more and more stable (or p2 no longer drops far below 1) as 
a is increased. 



This contrasts with the case of gA = 0.5, shown in 
recovers the situation of repulsively bound pairs [2, 



^ig.Ob). Near the static limit (a < 0.01), one 
1611 . which tunnel at a period T /2n ~ u/4J 2 



long compared with the Rabi oscillations 1/27, and which are stable in time (inset). Allowing for 
a finite mass of B, this pair breaks up, and the dynamics starts to resemble sinusoidal Rabi oscil- 
lations for a = 0.12. This becomes even more conclusive from the angle of the pair probability 
(inset), whose minimum value now significantly deviates from unity. Qualitatively, this feature is 
captured by the effective Bose-Hubbard model ©: At large enough a, the attractive on-site inter- 
action 8u < tends to cancel uq, which is in agreement with the Rabi-like oscillations observed 
here. 

What happens for increasing intra-species repulsion, in particular upon approaching the 
fermionization limit, e.g., at gA = 25 [Fig. He)]? In the case of a quasi-static effective barrier 
(a = 0.001), the dynamics resembles that known from fragmented-pair tunneling in a double 

n 

well [16]. By the Bose-Fermi duality, this may be regarded as two noninteracting fermions Rabi- 
tunneling independently in the lowest Na = 2 bands /3 = 0, 1: pr(?) = E/3 sin 2 (j^tj /Na- With 



increasing a, from our previous discussion we expect an effective attraction between these two 
fermions. In the intermediate regime (a = 0.01), this bears little effect and only leads to a renor- 
malization of the two "Rabi" frequencies However, for mass ratios as large as a = 0.12, 
the picture changes qualitatively: The two atoms tunnel only on a time scale about four times 
longer than the Rabi oscillations. This pattern closely resembles that of correlated pair tunneling. 
Indeed, a look into the two-body correlations [Fig.^c), inset] reveals that the pair probability P2 
stays remarkably close to unity, in marked contrast with the conventional fragmented- atom pair. 
In this light, it is enticing to think of this as the tunneling of an attractively bound pair of identical 
fermions. 

Let us account for this in a simplified model. Modifying the derivation of © for the case of 



noninteracting/ermzons, we find an effective multi-band Hubbard model II 1411 . For the special case 
of Na = 2 "fermions", this takes the form 

where Na bands /3 = 0, 1 contribute; note that we have discarded the on-site term £ s 
constant in our setup. The induced on-site interaction turns out to be negative, 8u < 0. For 
larger a, the induced interaction \8u\ 3> J^> shifts single-atom tunneling far off resonance. It 
is in this limit that we can understand the inter-species tunneling as that of an attractively bound 
fermion(ized) pair at a fairly large tunnel period T /2k ~ \8u\/AJ^'J^\ From this viewpoint, the 
two quasi-bound fermions reside in different bands, i.e, pertain to different pseudo-spins and thus 
are not constrained by Pauli's principle. This picture extends to arbitrary numbers of fermions 
Na > 2, which tend to tunnel in dynamically bound A^-atom clusters. 

In conclusion, we have investigated the tunneling of bosonic atoms through a second, localized, 
species. Remarkably, this can be well understood as an induced attraction between the mobile 
bosons. For the ground state, it may lead to strong intra-species correlations even in the absence 
of intra-species interactions. The inter-species dynamics, among other things, features tunneling 
of essentially an attractively bound fermionized pair. Extending the discussion of tunneling to 
nonclassical potential barriers opens up intriguing perspectives, such as creating an effective lattice 
made of localized atoms. This would allow for the study of disorder beyond the "quenched" static 
limit, to give but one example. 
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